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Based on the collective variables representation with the Hermitian form of Bogoliubov-Zubarev 
Hamiltonian the self-consistent oscillator model of a ground state and excited states of Bose liquid 
has been proposed. The new method of calculation of anharmonic terms in this Hamiltonian and 
its interpretation have been presented. The dispersion equation for a collective excitation in a 
superfluid ''He has been obtained in self-consistent way, where real and virtual processes of decay 
of a collective excitation were considered. The end point, determined by a threshold of collective 
excitation's decay on two rotons, of the dispersion curve has been obtained and it was shown that the 
dispersion curve strongly depends on property of its stability. An approach with a structure factor 
has been realized without using of any adaption parameters. Based on the oscillator model the new 
method of self-consistent calculation of a ground state energy and density of Bose condensate has 
been proposed. The model of suppression of Bose condensate has been presented. 

PACS numbers: 67.25.D-, 67.25. dt, 67.10. Ba 



I. INTRODUCTION 



Up to present time many various models have been proposed for microscopic description of a superfluid helium. 
Most of the models are based on two approaches. The first one is the theoretical-field formalism. Its characteristic 
feature is explicit application of Bose-condensate (BC) and introduction of the higher BC in some models [1-4]. The 
second one is the quantum-mechanical approach, where the Schrodinger equation is solved for a ground state and for 
lowest excited states of N interacting particles. This approach originates from the papers of Feynman and Cohen Q 
where connection between the structure factor and the dispersion curve of a collective excitation (CE) was obtained, 
however this connection is correct at small values of a wave vector only (/c — > 0). The quantum- mechanical approach 
does not use BC explicitly unlike the theoretical-field approach but it calculates BC knowing the wave function of a 
ground state. 

The quantum-mechanical approach has been developed essentially due to use of Bogoliubov-Zubarev formalism Q . 
In this representation Hamiltonian of a bosons' system is written in the terms pj^ and d/dp-^ where p-^ is a Fourier 
transform of fluctuation of density. As it has been shown in the papers ^, 's^ normal motion and superfluity motion 
can be separated in the iV-particle Schrodinger equation. The oscillations of Bose liquid is obtained in the harmonic 
approximation. Zero-point oscillations correspond to a ground state of liquid, and excited states - to any collective 
excitations: phonons, rotons, maxsons and so on. However the non-harmonic terms of the Hamiltonian play a essential 
role. The first correction to the harmonic approximation has been calculated in above-mentioned papers, the second 
correction has been calculated in the paper [9|. Unfortunately any small parameters aren't in this expansion and 
calculation of each new correction is accompanied by large mathematical difficulties. Contribution of the anharmonic 
terms in energy of a liquid has been evaluated by Brillouin-Wigner perturbation procedure in the papers TO, 11], the 
method of Green function has been used [l3 , the " shadow wave function" approach has been developed [13| . 

The dispersion curve of CE obtained by some authors [1, has satisfactory coincidence with the experimental 

spectrum if value of the wave vector is fc < 2.5^1^^. In the region k > 3A~^ the experimental spectrum reacher 
the "shelf, that is hybridization of CE with the two-roton level is observed fl^. The region k > 3A~^ has been 
investigated by the theoretical- field method in 15| , where it has been shown that the " shelf on the dispersion curve 
caused by decay of CE into two rotons and by the end of the spectrum. Calculation of the next corrections to the 
spectrum of CE doesn't give the end point of a dispersion curve. We have to note that the potential of interaction 
between atoms is known badly on small distances r < 2.5A. A ground state energy is very sensitive to interaction 
on small distances just. The situation is complicated by fact that He II is liquid however complete theory of liquid 
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doesn't exist till now. That's why authors proceed from gas representations or the models of quantum crystal |16l - ll8l |. 

Bogoliubov-Zubarev Hamiltonian is non-Hermitian because the transition from Cartesian coordinates to the collec- 
tive coordinates is non- unitary. This Hamiltonian has been used in the many aforesaid paper [^-0 [13, E^-H^ • But 
the non-hermicity must lead to violation of some theorems of quantum mechanics. However an operator of speed gj^ 
canonical conjugated with pj^ has been introduced in the paper of Sunakawa [IQ] , and his Hamiltonian represented 
by these terms is Hermitian. The Hermitian form of Bogoliubov-Zubarev Hamiltonian and the Jacobian of transition 
to collective coordinates have been obtained in the book [2^ . However the anharmonic terms of this Hamiltonian has 
not been calculated. 

In the Section superfluid motion and normal motion have been separated in the Schrodinger equation with the 
Hermitian form of Bogoliubov-Zubarev Hamiltonian. In the Section Hill the oscillator model of Bose liquid has been 
formulated in the Random Phase Approximation (RPA) . In the Section IIVI contribution of anharmonic terms of the 
Hamiltonian in the spectrum of CE has been calculated. It has been shown that this correction describes decay of CE 
and the dispersion curve has the end point kQ. In the Section |V] a ground state energy Eg has been calculated taking 
into account anharmonic terms in the Hamiltonian with help the oscillator model. The mechanism of suppression of 
BC has been described and density of BC at zero temperature has been calculated. 



II. THE EQUATION OF MOTION AND THE SPACE OF COLLECTIVE VARIABLES 

Let's consider N interacting Bose particles of mass m confined in a macroscopic volume V. Hamiltonian of the 
system has a form: 

^= E E (1) 

l<j<Ar l<'i<j<Ar 

where p is the operator of momentum of a particle, the operator <I>(|ri — rj|) is an energy of interaction of two particles. 
The waves functions must be symmetrical for any rearrangements of coordinates of any pairs from N particles and 
it describes oscillations of Bose liquid. The space of the collective variables is the suitable multitude of variables 
describing collective motion of a system [23, • Collective variables are Fourier transform of fluctuation of density 
An(r)=Ef=i^(r-r,)-iV/y: 

1 ^ 

= E cxp(-ikr) = pI- ipl, (2) 

.7 = 1 



where 

4 



N , N 



^^cos(kr), 4 = ^^sin(kr). (3) 



TV ^ VN . , 



And what's more the correlations have a place: 

Pk-P-k Pk = P-k' Pk = -P-k- (4) 

This means, that it is necessary to consider the values pj^ with indexes k from the half-space of their possible values 
only. 

The transition from the Cartesian coordinates (ri,...,rAr) to the variables pj^ is nonunitary because quantity of 
the Cartesian coordinates is 3iV but quantity of pj,-variables is infinity. Hence superfluous variables exist among 
collective variables. The transition to pj^-representation must be done with the Jacobian which equalizes a volume of 
the configuration space J dri... J dr^ = V'^ to a volume of the pjj.-space: 

^"^ = 11/ / ^'^py^ (5) 



where the prime at the symbol of multiplication means that variables k are taken from the half-space only. In a 
representation of the new variables ^ the Hamiltonian ([T|) has the form: 
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^(0) + ^ E ^(fci) (Pk,p-k, - l) > (6) 



and its name is Bogoliubov-Zubarev Hamiltonian. However this operator is non-Hermitian - the first term in the first 
brackets: p-^^ . This property of -ffg2 caused by the the nonmiitarian transition from Cartesian variables to 

collective variables p-^. 

In the paper 23] the Hermitian form of the Hamiltonian ([5]) has been obtained. Let the system is described by the 
wave function ijj which is normalized in pj,-representation: 



dpl / 

Let's introduce the wave functions normalized without the Jacobian J: 

= (8) 

Then Schrodinger equation is written as follows: 

H^^E^, H = J^/^H^r^J-^'\ (9) 

where the new Hamiltonian H must be Hermitian. Proceeding from this condition we can find the Jacobian J. Then 
we have: 



ki#o 

kiT^O ka^O ^1 kiT^O 



(10) 



where the Jacobian must be found from the equation: 



9 In J 1 '^'^'^"kik2 9 In J 

^ kiT^o 
Solution of this equation can be written as 

lnJ = lnC + ^— ^-f^=— ^ ... Y: Pq.-Pq„- (12) 

The constant C can be found from the condition ([S]). I our comprehension the equation ([SJ with the Hamiltonian 
(jlOp is motion equation in pj^-space. The equation (jlip is constraint equation in this space. In Cartesian coordinate 
we have a discrete system from N particles but in collective coordinates the system is regarded as continuum. The 
constraints equation removes the superfluous degrees of freedom. 

Normal motion and superfluid motion can be separated in the Schrodinger equation 7]. The wave function of a 
ground state describes zero-point oscillations of Bose liquid and has an exponential form e^. The whole wave function 
(with excited states) has a form: 

4^ = e^ip. (13) 

Then let's rewrite an energy of a liquid as 

E^Eo + E-Eo = Eo + Ee^t, (14) 
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where -Eq - is a ground state energy (energy of superfluid motion with wave function e^), E^xt - is the energy of 
excitation (energy of normal motion with wave function ip) . Then the equation ^ divides into the set of equations 
describing superfluid motion 



-E E 



dU dU 



e(ki,k2 



N 



dU 



E 

ki#o 



N 
2V 



dU dU 

= Eo + 



(15) 



E 



-£(fc) 



N 
2V 



Kfc) 



2V 



HO), 



and normal motion 



- E 

ki5^o 



9pk,9p-k, 

E\ - £(ki,k2) 



^_mj__di^ 
'^Pki ^P-ki 



dU dip 



The wave function -0 must be eigenfunction of a momentum operator 0, Q : 

d 



= Eext{k)ip 



■ E ^kiPk.-apj^ 



P-0 = P^. 



The wave function ■0 is eigenfunction of the momentum operator too: 

ji/2p = P0 =^ P0 = Pi/;. 

The proof of this fact is in Appendix [A] Then solution of the equations (fT5|) and (fT6|) must obey the condition: 



= 0-e' 



Pip = h{ki 



kn)ip. 



(16) 



(17) 



(18) 



(19) 



This means that a center of the system's mass rests and the function ip describes excited states of a system, which are 
characterized with conserving energy Eext{k) and with conserving momentum h{ki + . . . + k„). Thus the function 
ip describes state with n collective excitations. If an interaction of particles i^(fc) — then the wave function ip is a 
wave function of free bosons ipi,=a — l/Vv^. A proof of this statement is in Appendix iBl 
In the total case a solution of the equation ([T5|) has a form of correlation series expasion0, Q: 



k+q^^o 

c^ = E/w^kP-k + E E 

k#o k#o q^^o 



5(k,q) 

^Pq-k^kP-q + 



(20) 



Analogous expansion for the function ip of state with one CE (as Pp = hkip) has a form: 

k+q^a 

P-k + 



ip : 



E 



^(k,q) 

^^^q-kP-q 



(21) 



Unfortunately no small parameter exist in these expansions. And what's more, calculation of higher terms of the 
expansions is very difficult mathematically. In order to solve this problem we have to use some model considerations 
and approximations. 

We have to make a little remark hear. Under collective excitations we understand quants of collective motions of 
macroscopic group of particles, that is with motion of system as a whole. For example: phonon, plasmon, magnon. 
Under quasi-particle we understand a particle interacting with its environment or external field ("dressed" particle). 
The quasi-particles are characterized by effective mass and they are interacting with effective (screened) potential. 
The examples of quasi-particle are conduction electron, polaron. Cooper pair. 
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III. RANDOM PHASE APPROXIMATION 



A. The random phases as zeroth approximation 



In this section we will shortly formulate zeroth approximation to our problem - RPA or harmonic approximation. 
The Jacobian and the wave function of a ground state wc write in Gauss form: 



In J = InC - i ^ Pk^-k- 



(22) 



k^^o 



k^^o 



(23) 



Substituting these expressions in the equation (ITSI) and neglecting by powers of pjj. higher than second power we have: 



-4£(fc)/2(fc) + i£(fc) + ^Kfc) 



E 

Eo - ^KO) + E 
k^^o 



PkP-k 



2e{k)f{k) + ie(fc) + ^u{k) 



From this equation we can obtain unknown function /: 



f{k) 



2N v{k) 



and corresponding energy: 



k#o k/o 



V e(fc)' 



1 N 
-e(k) H iy(k) 

2 ^ ' 2V 



7V2 

—iy(0). 

2V ^ ' 



It is necessary to notice that the solutions 



correspond to the shortened Hamiltonian: 



a2 



lpk,p-k. 



2V 



where, unlike the complete Hamiltonian ()10p . the harmonic terms is kept only. 
A structure factor is the most important characteristic of liquids [2^ : 



(24) 



(25) 



(26) 



(27) 



(28) 



S{k) ^ (Pk/^-k)' 

where the mean value () is calculated with a ground state. This function can be obtained by using the virial theorem: 

6E / SH \ 1 1 



Siy{k) \ 5v{k) 



S{k) 



1 



2N fjk) 
V e(k) 



4/(fc) 



(29) 



For atoms of helium the potential of interaction is very like to Lennard- Jones potential in the region r > 2.5A. 
In the remaining region the interaction is known badly and various adaption functions have be used here. In order 
to get over this difficulty the approach with structure factor has been developed in some papers where the potential 
of interaction is unknown function and the structure factor is taken from experiment data. This means, that the 
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perturbation theory must be constructed so we must obtain the same structure factor in each approximation. On the 
present step of calculation the potential of interaction is obtained from (1291) in a form: 



i^(k) = —e(k) 



1 



52(fc) 



1 



In order to obtain dispersion curve of a CE we assume that ip = Substituting its and the function 

equation (fT6|) we have: 



E. 



RPA 



e{k) / .,.„ 2N ,,, 



(30) 
in the 

(31) 



This expression is known as Feynman formula or Bogoliubov specter. The formula (j3ip is a starting-point in order to 
obtain the dispersion curve of a CE which is consistent with the experimental specter. 



B. Formulation of the oscillator model 



The Hamiltonian (|27|) is sum of separate terms where each of them is characterized by own wave vector k. This 
fact means that infinity number of independent motions can be in a system, and every of them is characterized by 
wave vector k and some energy £,(k). For the each independent motion Schrodinger equation is written as follows: 



9V 



2(TO/fc2) dp-^dp -^ 



-e{k) 



N 
2V 



This equation is similar to an equation of a harmonic oscillator, where 



2 1,2 



2m 



1 



N 



1 



^(fc) is energy of a oscillator with the wave vector k: 

m = hL^ik) Q 

and the ground state energy of a system is 

Eo = Y^ Ihcik) - J2 

k#o k#o 
The wave function of a ground state has the form: 



7eik) + —i^ik) = -Afw 



g(fc) (]_ 
S{k) \2 



1 N 



, Pkp-k\ /(fc)pl,p 1, 

expl ^'=^ 



Muj 



where I is a oscillator length. The structure factor can be found from the virial theorem for an oscillator: 



(PkP-k) 



s 



nuj{k) ' 



2ln) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



and the state n — corresponds to a ground state. A structure factor S{k) is dispersion of a oscillator with a wave 
vector k. 

On the assumption of the aforesaid we can formulate the oscillator model of Bose liquid. A quantum liquid represents 
a totality of harmonic oscillators. Each oscillator is a vibration mode of density of liquid. A ground state energy of 
a system is a sum of energies of ground states of these oscillators. Availability in a system n collective excitations 
corresponds to a state with n excited oscillators. The oscillator model is a model of continuum medium and it is 
incorrect on distances between neighboring atoms. The calculation of short-range correlations is very difficult because 
all anharmonic terms in the Hamiltonian ([T0| and in the Jacobian ([T2|) must be calculated. Moreover if we shall 
try to calculate them with help ordinary perturbation theory then we shall have divergent integrals. This difficulty 
appears due contribution of above-mentioned superfluous degrees of freedom. In the next sections we shall formulate 
the method of calculation of the anharmonic terms remaining in limits of the oscillator model. 
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Figure 1: The spectrum of CE calculated with the Feynman formula (solid line) and measured in experiment (dotted line). 
Decay of excitation on two rotons is possible in the region k > where E^J^^{k) > 2A^^^. The experimental curve has 

a "shelf at fc > 2.5A^^ where decay of CE into two rotons is observed. 



IV. THE SPECTRUM OF COLLECTIVE EXCITATION 



A. Decay of collective excitation 

The wave function of a excited state ip — leads to the well known Feynman formula (1311) . We can see on Fig[T] 
that the dispersion curve in RPA coincides with the experimental dispersion curve in the limit fc — > only. Its basic 
contrast to the real spectrum is that that Feynman formula doesn't take into account decay of CE and the end of a 
dispersion curve accordingly. 

Let's consider the wave function of CE ((2T|) with momentum Kk where L(k, q) is the unknown function. In our 
interpretation the second term in (|21|) describes a process of decay of CE with momentum Kk and energy Eextik) into 
two excitations with momentums ?i(k — q) and ?iq and with energies Eext{\i^ — q|) and Ef,xt{q) accordingly. Decay of 
CE can occur if energy is conserved: 

Eextik) = £;exi(|k - q|) + E,xt{q). (38) 

Such situation corresponds to the case (A) on Figl2] This process causes damping of a CE and the end of the 
dispersion curve [1^. In those parts of the energy spectrum where the equality (pS)) isn't true, CE is stable. In this 
case the term in (|2ip with the multiplier P(^_\i^P-<\ describes processes of virtual decay and creation of CE. This 
causes trivial renormalization of the energy spectrum of CE. This situation corresponds to the case (B) on Figl2] 
Higher approximations in ([21]) describe decays of CE on three and more excitations. We suppose that such processes 
is improbable and we shall not consider them in this paper. 

Substituting the wave function (|21l) in the equation (|16p we can obtain the function 



w. . £(k,q)g-i(g) 

E^xtik) - Ef^,^{\k - q|) - Ef^,\q) ' ^''^ 

where E^J^^{q) is spectrum of CE in RPA - the formula (PT|) . And besides we can obtain the dispersion law of CE 
in a form: 

E^xtik) = + 1 E - q) 

q#o 

- FRPA(h\ - — £(k, q)£(k-q,q)5'-l(g) 

N^^Eff,Wk-c\) + E^^P,\q)-E.xt{ky ^^"^ 

This equation is an integral equation for unknown function Eext{k). The first term corresponds to the term p_\^ m 
the wave function (j21l) and it is Bogoliubov-Feynman spectrum. The second term has a singularity in those points q 
where the equality is executed: 

E^xtik) = E^j;^i\k q|) + E^j;^iq). (41) 
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Figure 2: The processes described by the second term in the wave function (|21[) . The process (A) is a real decay, the process 
(B) is a virtual decay of CE with momentum ?ik into two excitations with momentums ?i(k — q) and ftq. 



This equation likes the expression (|38|) and it is condition of decay of CE. Inexactitude of this condition is in that 
that a real CE Eext{k) decays on excitations with a spectrum calculated in RPA: E^^^[q). However in the region 
of large q the theoretical dispersion law of CE differs from the real law Eext{q) essentially. We shall get over this 
difficulty later. For the present we shall work with the condition (|4T|) . 

CE is stable at those k where the equation (|^T|l hasn't solutions. Damping of CE appears in a threshold of decay 
k = k(j. At the point k = kn a solution of the equation ([^T|) appears at first, and the right sides of the equations ([55]) 
and ([4T|) have extremum [2^ as a function of q. We can see on Fig[T]that the extremum for the function Ef^^{q) is 
the roton minimum. In its neighborhood the spectrum of collective excitations has a view: 

i?f,^^(g) = A«^^ + |(<Z-M^ (42) 

(Ik - q|) = A«^^ + ^(|k - q| - fco)^ (43) 

where A^^'^ = is the roton gap in RPA (in degrees of Kelvin), ^ = 0.09m is effective mass of a roton in 

RPA, /co is momentum corresponding to the roton gap (in angstroms). We can suppose that fco = fc^^^ = 1.93^-1 
in all approximations. The condition (|4ip is executed when Eext{k) = 2A^^"*. Then the underintegral function in 
the spectrum (|40|) at k — kc has a singularity in the point q = qo- This means decay of CE with momentum hk(j 
into two rotons with momentums hko each. Hence we have the connection k(j — 2ko cos 9c where 29c is a angle of 
recession of two rotons. 

Let's evaluate the integral in (l40t in the point k = kQ. Let's consider asymptotics of the underintegral function at 
q — >■ oo where S{q) — >■ 1. If a vector k is directed along the axe Oz then kq — kqcosO. Let's use the expressions 
and (^5)) in the denominator of as they determine the singularity at g = fco, and we have |k — q| « q. Then the 
integral in (j40p is written in a form: 

,2 s 2 



V f fi \ ^ f {qk COS 9yq^dq sin 9d9dip 



N(2^f\2mJ V 2ARPA~Ecxt + ^{q-koY 
qk cos 9q'^ q^ dq sin OdOd^p 
2ARPA_Ecxt + ^{q-k^Y^ 

V f ri^V 2^^ f q^dq 



(44) 



/ ^ ^ ("45) 

N{2TT)^\2mJ 3 J 2ARPA-Ecxt + ^iq^kor' 

In the second integral in integration by 9 gives zero. Taking into consideration Eext{k(j) = 2A^^^ we can rewrite 
(|45| as 

2 V fh'\'2 f, r l' . ^ 

T2 I 7Z hrv?."'1 ~^ °° (46) 



iV(27r)2 \2mJ 3 (q - ko) 



The underintegral function has a pole in the point /cg that results in some complex addendum in an energy of CE: 
Eext — i/t. This means damping and limited life time ^ r of CE. Moreover we can see that ultraviolet divergence takes 
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place. Hence the second term in the formula (|40)) for the spectrum of CE aspires to infinity: Ei,xt{k) = E^^^{k) +00. 
The nature of this divergence lies in the following. 

The constraint equation removes the superfluous degrees of freedom appeared at nonunitarian transition to 
collective coordinates. These degrees of freedom are connected with continuum representation of primarily discreet 
medium. Each superfluous degree of freedom gives contribution to energy. Neglect of all anharmonic terms in the 
Jacobian (fT2|) and in the Hamiltonian (fTO|) doesn't give possibility to delete this nonphysical energy. Contribution of 
the superfluous degrees of freedom results in infinity in the integral (j40p . 

The method of removal of the divergence is cut-off of the integral (|40p to some unknown value of a wave vector 
Q ~ Qm- By analogy with [l5l. [26| let's introduce new variables q'^ q'p in accordance with the definition: 

qx = (fco sin 9c + q'p) cos </? qy = (/cq sin Oc + q'p) sin ip 

~ h) cos 9c + q'^_ (fiq = q^dq sin 9d9dip = (ko sin 9^ + q'p)dq'pdq'^dip. (47) 

In a neighborhood of the threshold of decay k = k(j (where Eext 2A^^^) the imderintegral function in PO)) 
has a pole in the point [q = ko, 9 = 9c]- We have \q'J <C fco and \q'p\ <C fco in a neighborhood of the singularity 
g — )■ ko,9 ^ 9c- Then we can write expansions near the singularity: 

d^q ~ ko sin 9cdq'pdq'^d'-p, q ko + q'^ cos 9c + q'p sin 9c 

|k - q| w /jo + q'pSindc - q'z cos9c, kq « kkocos9c + 0{q'p,q'J 

(k-q)q^klcos29c + Oiq'p,q',)- (48) 
In the limit q ^ k^ energies of products of decay have a view: 

Sf,^ (g) = A«^^ + ^^(qi COS0, + q'p sin9c f, 

(|k - q|) = A^^-^ + ^i<l>^(^c - q'z cos9cr- (49) 

Let's consider the spectrum of CE (|40p in some neighborhood of the threshold of decay fc — /cq. Decomposing the 
underintegral function near the pole q ^ ko, 9 ^ 9c using the expansions and (02) we have:: 

2 ^ e(k,q)£(k-q,q)5-l(g) ^^^^ 



N ^„ i?«^^(|k - q|) + E^J^^iq) - Ecxdk) 



V (n^kiy n 



2 

^"^(fco) cos 29c cos 9c sin( 



N (27r)2 \2m J 

f dq'^dq'p 

' ^{2ARPA^Ecxt)+[{q'p}^sm^0c + {q',ycos^9c]' ^ ' 



However the integrals ([501) and (jSTt are not equal because transition to the limit q ^ ko, 9 ^ 9c has been done 
in the underintegral function. Integration by all g-space results to infinity. Then we must introduce some cut-off 
parameter p so as the integral ([5T|) is finite on the one hand, and the integrals ([50)) and (|5ip are approximately equal 
in a neighborhood of the point q = ko: ko — p < q < ko + p on the other hand. 
Now let's introduce the polar coordinates: 

q'p sin 9c = p cos ip 
q'^ cos 9c = p sin ip 

pdpdjp 

dqpdq^ = n ■ n ' 

cos 9c sm 9c 

Then the integral (ISTI) is reduced to the form: 

_ 2 y Mcc^ r 7^ — (53) 

^ \^f{ko}icos20c^^^ \ {2APP^ - Ecxt)fi/n^ 
~ N 2T:fi^S{ko) 



{2ARPA - Ecxt)t-i/h' + p 



10 



We can see that the final result depends on the unknown parameter p which determines a integration domain near 
the point q = kQ. Then p must satisfy the requirements: 

1. The point k(j is the threshold of decay: Eext{kQ) — 2tJ^^^ 

2. On FigUwe can see that Feynman formula E^J^^{k) — e{k)/S{k) in the limit fc — >■ gives a sound mode chk 
with a correct sound velocity c. Hence the total spectrum must have the same asymptotic: 

E,,t{k^O)=E^P^ik^O) = chk. (54) 

The second term of the expression (j^D]) in the limit k ^ k(j is the integral ([51]). We must join these two 
asymptotics in a single expression which is a sought dispersion curve. Hence the parameter p is function of 
variable k: p ^ p{k)- 

3. The parameter p must ensure approximate equality of the integrals (1501) and (I5ip in a neighborhood of the point 
q = kg. We have decomposed the underintegral function in (I50|) near the point g = fco in power series of p/kg: 

(J « fco + 9z cos^c + q'p sinOc « fco (1 + p/ko) , 
|k- q| sa fco + gpSin6'c - q'^ coaOc ~ fco {I + p/ko) , 

kasmOc + q'^ kosin0c(l+ , ^„ ) ^ , . < 1- (55) 
\ ko sm y ko sm t^c 



Hence p(fc) must satisfy the inequality ([51 
4. The spectrum of CE (1401) must be determined with the parameter p{k) by self-consistent way. 
In order to meet the requirements we must assume that 



(2A«^^ - Ee.t)t^/h^ _ 2xA«^^ - E, 



ext 



i2A^P^ - E,,th/h' + p^ 2xA«^^ 

then 



(56) 



2,, , _ M ge.t(fc) [2A^^^-i^,.t(fc)] 

^ - ¥ 2xA^P^ -E^^,{k) ' ^^^^ 



where the parameter x is determined by the condition ([T]) Eext{k(j) = 2A^^"^: 



As a result the dispersion equation for CE (|40p has a form: 

2xA«^^ - E,,tik) 



E,xt{k)=EgP'^{k)+ak\n 
where 



2xA-R^^ 



(58) 



(59) 



a=— ^ , cost'c = 7rr^. 60) 

We can see that the parameter p is determined by the spectrum E^xt (k). On the other hand, p determines this spectrum. 
Eext{k) doesn't contain p in an explicit form. Hence Eext{k) is determined with the parameter p by self-consistent 
way. 

In the long- wave asymptotic E^^{k — > 0) = cUk + Ak^ , A — const > 0. Then it is easily to show that 

Eextik 0) ^ chk - ^^^RPA k'+Ae-.... (61) 

This means that the spectrum E^xt is joined with Feynman formula at small k and describes the sound mode. Thus 
the condition ^ is satisfied. However we can see that even powers of k exist in the expansion (pT|) . But this 
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E„,(K) 



2A 



(1) (4) (7) 



2A'^^'^ (RPA)/ // 















ma") 



Figure 3: The spectrums of collective excitations. The line (RPA) is spectrum corresponding to Feynman formula (|3ip . The 
lines 1-7 are dispersion curves obtained as numerical solution of the equation (|59p with parameters fcQ equaling 2.9yl~^, 3.0A~^, 
3.1y4~^, 3.2j4~^, 3.3A~^, 3.4A~^, 3.5yl~^ accordingly. The dotted line is experimental spectrum of CE ^] at low temperature. 
The dispersion curves 1-7 finish in the points k = /cq on the line 2A^^'*. 



P(A") 




/C(A-') 



Figure 4: The parameters of cutt-off p(fc) for the integral (|53|) . The curves p(fc) correspond to the curves 1-7 on Fig. ((3]). A 
domain of the functions p(k) is < fc < 



doesn't contradict to isotropy of liquid because the dispersion equation (j59p depends on modulus of wave vector only 

Eext = ^^ea;t(|k|). 

We can see that the dispersion curve ([55]) depends on the parameter k(j (threshold of decay) strongly. This is free 
parameter being in the limits < k(j < 2fco, therefore we have a family of curves shown in Fig|31 Corresponding 

cut-off parameters p which joins the long-wave asymptotic /c — >■ with the asymptotic k — ^ k(j is shown on FigU 

A domain of the function p{k) is < fc < kfj because at fc > kfj we have p{k)'^ < 0. This means that the dispersion 
curve (j59p can not be continued to the region fc > k(j (though the function Eext{k) exists formally). Hence the point 
k(j is end point of the spectrum. For the typical value kQ = 3.2^4^^ we have: max(p/fco) = 0.28/1.93 = 0.15 <C 1 and 
max{p/ko/ sin^c) = 0.27 < 1. This means that the condition ([3|) is satisfied. 

The dispersion curves 1-7 on Fig 13] are more close to the experimental curve than E^J^^{k), however they have 
an essential defect. CE with the dispersion curve E^xtik) decays into excitations with spectrum calculated in RPA: 
^ext^il) E^^^{\\i — q|). This situation corresponds to conservation of energy in the form In reality 

the law of conservation of energy is the form ([38]) where energy of decaying CE and energies of decay products are 
determined with the real spectrum. We can see on Figj3] that the dispersion curve Eext{k) differs from the curve 
E^J^^{k) essentially in the region of big fc. So, the new value of a roton gap is A < A^^"^. This means that energy 
is not conserved at decay of CE into two rotons. In other words the dispersion curve Eext{k) is not self-consistent. 
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Figure 5: The spectrums of collective excitations. The line (RPA) is spectrum corresponding to Feynman formula (|3ip . The 
lines 1-7 are dispersion curves obtained by numerical solution of the set of equations (|62|) and (|64p with parameters kQ = 2.9A~^ , 
3.04"\ 3.1A"\ 3.2^"\ 3.3^"\ 3.4A"\ 3.5yl"^ accordingly. The dotted line is the experimental spectrum of CE [2^. The 
dispersion curves 1-7 are ended in the points fc = fcQ on the lines 2Ai_7. 



B. The self-consistent form of the dispersion curve Eext(k) 



The dispersion equation E^xt (k) is self-consistent if an energy of the end point of the spectrum is equal to the doubled 
Toton gap of this spectrum: EgxtikQ) — 2A. Let the new spectrum has a form 



Eext{k) 



2A-E^xt{k) 
2A 



where the parameter A must be such that the condition of self-consistency is executed: 



EextikQ) 
Eext{ko) 



2A 

: A 



The parameter A is the new roton gap in the point /cq. The the set of equations (|63p in a expanded form is 



2A^E^f,^{kc)+akc\n(l 



A^E 



RPA 



(fco) -I- afco In 1 + 



2A 



(62) 



(63) 



(64) 



The set of equations ((62|) and ((64)) determines a dispersion curve of CE. The results of numerical solution of these 
equations with kQ = 2.9A^^ — 3.5A~^ are shown on FiglS] The curves 1-7 on FiglS] correspond to the curves 1- 
7 on FiglH As in the previous case the spectrum depends on its end point k(j strongly though sound velocity 
he = limfc_j.o is the same in the all cases. We can see a weak pinning of the curve Eext{k) to the line 2 A. 

This means hybridization of the dispersion curve of CE with a two-roton level [27]. On the rest of regions an energy 
Eext{k) is determined by a state with one CE (phonon or roton) completely. We have the family of curves Eext and 
we cannot obtain the point k ^ on t his step of calculations. In order to do this we must generalize the oscillator model 
formulated in the Subsection IIII Bl 



V. THE GROUND STATE AND BOSE CONDENSATE 



A. Formalism of effective mass in the oscillator model. 



In this section our problem is calculation of a ground state energy taking into consideration the anharmonic 
correction in the equation p^ . If we calculate this corrections then we have a divergence Eq = E^^^ -I- oo as in the 
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Figure 6: The effective masses of quasi-particles in units of mass of "^He: m/m. The curve (RPA) corresponds to a "naked" 
particle. The numbers 1-7 correspond to curves in FigO The functions m{k) have a sense in the range k < kfj only. 



Section ITVl In order to overcome this difficulty let's use the oscillator model formulated in the Section Hill where a 
superfluid liquid is a totality of independent harmonic oscillators with frequencies E^^^{k) = e{k)/S{k). In higher 
approximations we have another spectrums of CE Ee^itik): the set of equations (|62p and (|64p. However this spectrum 
can be written by analogy with RPA as 



Eext{k) 



m 

s{ky 



where 



e{k)^ 



2 1,2 



2m{k)' 



(65) 



Due to contribution of anharmonic terms of the Hamiltonian (jlOp mass of a particle m is renormalized to mass of a 
quasi-particle fh(k). Continuing the analogy with RPA we can write the effective Hamiltonian having a form as in 
RPA: 




^eff = ^(^i) 



(66) 



ki#o 



and the corresponding Jacobian: 



In J = In 



\ k#o / k^o 



(67) 



where the normalization requirement ([5]) has been used. 

We must obtain the same structure factor S{k) in any approximation because this function is given by experiment. 
Therefore we must introduce the effective interaction i>(fc) as follows: 



S{k) 



4/(fc) 



1 I 2N lyjk) 
V e(fc) 



v{k) = e{k) 



V 
2N 



S{kf 



2N yjk) 
V e(k) 



(68) 



The effective masses and interactions are shown in FiglB]and FiglT] It is necessary to notice that fn{k — > 0) = m and 

m = z.(o). 

As a result we have the totality of harmonic oscillators again. These oscillators have the same dispersions {p]^P_\^ 
as in RPA. However they have another frequencies E(,xt{k) instead of the old E^J^^{k). Difference between new 
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Figure 7: The potential of interaction between particles (RPA) and between quasi-particles (1-7). 



frequencies and the initial frequencies depends on the aforesaid anharmonicities. A motion equation for each oscillator 
can be written as the equation (1321) but with an effective mass and an effective interaction: 



2M{k) 5Pk"9P-^k 
A ground state energy can be rewritten as follows 

1 i{k) 



N 
2V 



2S{k) 



E 

k^o 



—V(k) 
2V ^ ' 



7V2 

— ?(0). 



(69) 



(70) 



It is necessary to notice that we must do transition from summation on k to integration as follows: 



k^^o 



V 



(2^)3 



k^dk. 



(71) 



The results of calculation of a ground state energy (per one atom) for particles' system (RPA) and for quasi-particles' 
system (the curves 1-7) are: 



RPA 



/N = -Vi.7%°K, e''^'^ /N = -^.7TK, E^^^' /N = -6M°K 



(2)/w_ 



E^^^ /N = -A.62°K E^"^^ /N ^ ~2m°K E^^^ /N = -0.59°K 



Ef^ /N = 1.43°i^ E}^'' /N ^ 3.43°i^ 



(7) , 



(72) 



We can see that existence of the end point of a dispersion curve causes some increase of a ground state energy. Since 
kQ is unknown then we cannot select a value of this energy. The energy Eq is minimal at kQ < 2.9yl^^. But the 
important condition for oscillator frequencies huj = E^xt exists yet. 



If velocity of a CE 



dp 



is more than a sound velocity in a system then CE radiates a phonon 26]. The 



regions of a dispersion curve, where such is possible, is unstable. We have formulated the oscillator model of Bose 
hquid where its ground state energy is a sum of energies of ground states of oscillators representing the liquid. In 
other words the ground state is an infinity number of virtual collective excitations. Then in order to determine a 
ground state as a state with minimal energy a dispersion curve E^xtik) must he stable. It means that for all k the 
inequality must be executed: 



dEextjk) j.^^ dEextjk) 
dk fc-i-o dk 



ch. 



(73) 



If this condition is not satisfied then the ground state falls apart. The plots of CE's velocity is shown in FiglS] The 
curves 4-7 satisfy the criterium ([73|) . Hence possible end points present in the interval k(j = 3.2A~^ 2ko only. The 
spectrum 4 corresponds to the minimal value k(j — 3.2^^^. On the contrary, the curve (RPA) is unstable almost 
completely. 
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Figure 8: The velocities of collective excitations Vext{k) = ^^na'k'''^ P®'' units of a sound velocity c. Vext{k) < means that the 
velocity and the momentum k have contrary directions. The curve (RPA) corresponds to Feynman spectrum. The curves 1-7 
correspond to the curves on FigO The condition (|73l) is satisfied for the curves where \vext{k)\/c < 1 for all k. The spectrums 
4-7 are stable about radiation of a phonon. The spectrum 4 (with k(j = 3.2A~^ , it is marked by bold line) is on the limit of 
stability. 




Figure 9: The dispersion curve of CE obtained with the optimal parameter k(j = 3.2 A ^. As result we have: fj, — 0.19m (the 
new roton mass), A = 8.S°K, ko — 1.95yl~^. The dotted line is experimental spectrum of CE [2^ at low temperature. In the 
point k = kQ on the line 2A the dispersion curve ends. 



We can see in (TfS]) that increase of kfj causes increase of Eq: Eq < Eq ' < Eq < Eq ' . Hence the energy Eq is 
minimal for all allowed energies Eq . This means that 

Eo/N = E^Q^^/N = -2.60°^. (74) 

The experimental value is Eq/N — ~7.12°K. It is necessary to notice that short-range correlations make essential 
contributions in E^. These correlations are calculated badly by the collective variables method. So, in [l^ it has 
been purposed the hybrid approach where transition to collective variables was made not completely and part of 
the coordinates was kept as Cartesian coordinates. Proceeding from aforesaid we can consider the result (17^ as 
satisfactory. 

The final spectrum of a collective excitation E^xtik) is shown in FiglH] The dispersion curve has been obtained 
as a numerical solution of the set of equations (p^ and with the optimal parameter k(j = 3.2 A~^. We can see 
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a weak pinning of the dispersion curve to the hne 2A that means decay of CE into two rotons with energy A and 
momentum fco in each. The dispersion curve coincides with the experimental curve weU. 



B. Suppression of Bose condensate. 



In Bose system at sufficiently low temperature the occupation number Nq of a state with minimal energy (momentum 
of particles is p = 0) is macroscopic number Nq < N. For ideal gas density of a condensate and temperature of the 
transition are 



No 
N 



1- ^ 



3/2 



To 



27rh^ ( N/V 
2.612 



2/3 



(75) 



At temperature T = all atoms are in BC: Nq = N . If particles interact then BC is suppressed. BC is suppressed 
strongly in Hell: Nq/N — 0.07. In the paper 3,j the model was proposed where BC is suppressed due to formation of 
pairs of boson. In this section we shall consider suppression of BC from the standpoint of collective motions of liquid 
using the oscillator model. 

The equations ([5^ or describe harmonic oscillations of Bose liquid. According to the uncertainty principle 
the energy of a oscillator in a ground state can't be zero \/2Tiuj ^ 0. The energy is divisible by kinetic energy and 
potential energy in equal parts (T) = ([/) = Afcj^(/9]^p_]^)/2 — £_/2 according to the virial theorem. This zero point 
energy is the energy of fluctuations of a quantum liquid. Hence collective motions is in a ground state and no all 
particles can be in BC. That's BC is suppressed by the dynamical fluctuations (29j . 

Using the equation ([M)) let's write the energy of a ground state in the a form: 



-4/2(fc)pkP_k - 2/(fc) + \pkP-k - I 



N^ ^ N ^ 
'^'"o + ^"kiPhP-k - 1) = ^0 



(76) 



After ground state averaging of this expression according to the rule P]^p_\^ ^ (P'kP-'k)^ taking into consideration 
the expression for the potential energy of a system 



(77) 



and the evident equality Eq ~ (T) + ($), we can write the average kinetic energy of a system as 



(r)-E^(^) 



'^f{k) + -\ (pkP_k) - 2/(fc) - \ 



(78) 



Since e(fc) is kinetic energy of a particle then the expression in the square brackets can be understood as occupation 
numbers of above-condensate particles N-^^^. Using the connection ()29|) we can write the kinetic energy and the 
corresponding occupation number via the function f{k) or via the structure factor S{k): 



w = E 



2 1,2 



2m 



(4/(fc) + 1 



\21 



-16/(fc) 



N 



k^o 



k^o 

{4f{k) + 1)2 (Sik) - 1)2 



-16/(fc) 



Then the density of BC is 



iVo = TV - E = ^ " 



V 



(27r)3 



AS{k) ■ 

(^-1)^ 
45* 



Ank dk. 



(79) 
(80) 

(81) 
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In the case of non-interacting particles we have: /(fc) = —1/4 and S{k) = 1, hence N-^^ = 0, Nq = N. In the case of 



weakly non-ideal gas with the interaction t^j^ = ^-^^a if a scattering length -^^r^ ^ 1 the density of BC is 



N 3 



(82) 



that coincides with the result of Bogoliubov ^26] . 

The expression for density of BC (|81l) is correct for systems with weak interaction only. As it has been shown in 
[20I [2l| for a total case we must use the expression 



iVo 
N 



F{R ^ 00) — exp 



V 



N^Ank^dk 



(83) 



where i^(r|r') is a one-particle density matrix, R — |r' — r|. The expression (|8T|) is the two first terms of an expansion 
of the exponent ([83|) . In RPA the density of BC is 



No 
N 



0.263. 



(84) 



BPA 



In the higher approximations we have the effective Hamiltonian (j66| . In order to calculate density of BC with the 
effective Hamiltonian H^q (|66p we are going to proceed from the following reason. Let's write the energy of a ground 
state of Bose system in a form: 



(5(fc) - 1)= 



45(fc) 



^"o + l^^'^kiSik) 



k#o 



2V 



Eo 



(85) 



We must mark out the kinetics energy of a system of particles. All anharmonic terms in the Hamiltonian if (|1QD or in 

.0 

the Hamiltonian ^^32 ® terms of a operator of kinetic energy J^i 



Ll 

<j<N 2m 



. This means that renormalization 



of mass m ^ m and interaction i/{k) — > i/{k) occurs due to the contribution of kinetic energy but not the interaction 
as in usual perturbation theory. Hence the potential energy of a system J2i<i<j<N ~ (EH) does't change at 
transition from H to H^q. In order to obtain the new kinetic energy of particles let's rewrite (j85p in the identical 
form: 



2V 



iSjk) - 1) 
45(fc) 



+ ^{^k-'^k)[S{k)^l 



(86) 



Proceeding from aforesaid the first term is a new kinetic energy of particles. In order to obtain a momentum 
distribution iVj^ of particles let's mark out kinetic energy of a particle e{k) in this expression. Then we have 



{T) = Y.e{k) 



m 

e(k) 



iSjk) - 1) 
4S{k) 



2V e{k) 



S{k) 1] 



(87) 



e{k) 



{S{k) - 1)2 



45(fc) 



2V e{k) 



S{k) 1] 



Substituting ((55)) in ((55)) we obtain the densities of BC for the dispersion curves 1-7 in FigjSJ 



iVo^^V^ = 0.149, iV^^ViV = 0.104, iVo^'^ViV = 0.071 
N^^Vn=0.048 



r(2 



(3), 



N^^^ /N = 0.031 



N^^^ /N = 0.020 



N^^'^ /N = 0.013 



(89) 



As it has been said before, the curve 4 with /cq = 3.2A~^ realizes in practice only. This means that portion of BC 
in superfluid helium at zeroth temperature is Nq/N = 0.048. Density of BC measured in experiment by the method 
of deep inelastic scattering of neutrons [s^ and by the method of quantum evaporation [sij is 0.07. Hence the result 
is satisfactory fully. 
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VI. CONCLUSION 



As a result of self-consistent solution of Schrodinger equation for the state of Bose liquid with one CE the dispersion 
law Eextik) (the set of equations and (|64l) ') has been obtained. It has been shown that for the Hermitian form 
of Bogoliubov-Zubarev Hamiltonian the second term in the Bijl-Dingle-Jastrow expansion considers real and virtual 
processes of decay of CE. Higher terms in this expansion considers decays into three and more excitations that is very 
improbable processes. Decay of CE on two rotons causes existence of the end point fcc of & dispersion curve. Our 
main result is: the phonon-roton dispersion curve E^xtik) determined by both interaction between bosons v{k) and 
the end point kQ . That is the dispersion curve strongly depends on property of its stability. 

In order to find fee and the energy of ground state we have been formulated the oscillator model of Bose liquid. 
According to this model the basic B-Z Hamiltonian can by rewritten as Hamiltonian in the harmonic approximation 
with renormalized kinetic and potential energies. As a result we have a totality of harmonic oscillators again with 
the same dispersions >5'(fc) but with another frequencies already. Energy of a ground state depends on these 
frequencies, hence it depends on the end point Uq. Proceeding from the condition of minimum of energy and from 
the condition of stability of ground state we can obtain the unknown end point = 3.2A^^ of the spectrum. The 
obtained spectrum of CE E^xt is a function connecting energy of a excitation with the structure factor S{k) of liquid 
which is taken from the experiment. Thus we did not use any model potentials of interaction. It means that all our 
calculations have been done without any adaption parameters. 

A mechanism of suppression of one-particle BC has been described at temperature of absolute zero. As a consequence 
of the uncertainty principal the dynamical fluctuations exist in a ground state and give nonzero kinetic energy. These 
fluctuations suppresses BC. It has been shown that such fluctuations are absent for ideal gas. It is necessary to notice 
that the analogy dynamical quantum fluctuations in other systems can cause quantum phase transitions (3^ . At 
nonzero temperatures BC is suppressed by both dynamical fluctuations and kinematic (thermal) fluctuations. 



Appendix A: The operator of momentum 



Let's consider the expression: 



/2. 



ki5^o 



d 



k^iokiTto 



? E E E ^^^P\<i^ 



d 



k/0 q#oki#o 



gPk^q-k^-q 



+ . . . + P^/; = - i ^ [?ikpkP_k - ^kpkP_k] + J E E [^l^/'kPq-k/'-q 
k^^^o k#o q#o 

+n(q - k)pk/Oq_k'°-q + H-^)P]^P(l-^kP-^\ - • . • + = + P-0 =^ P-0 PV5. 



(Al) 



That is if the function -0 is eigenfunction of the operator of momentum P then -0 is eigenfunction of the operator of 
momentum too. 



Appendix B: The wave function of a system of free bosons 

If to exclude interaction between particles: ^{k) = 0, the wave function of a system -0 = 0J~^/^ must be a wave 
function of free bosons. So, we can write for a ground state: 

= Be^J-^l'^ = const =^ [/ = i In J, (Bl) 

where i? is a normalization constant. Let's verify it. In a ground state of a system of nonintcracting bosons E^ = 0. 
Then proceeding from the equation (jisp we can write: 

d'^U dU dU 1 dlnJ 1 

(B2) 
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klVN 



dU 



dU dU 



= 0. 



Now let's rewrite the third term of the equation (jB2p in the form 



1 ain J 



4 9pk, 



^/'-kl 



J_ V — 
TTff ^ T2"^k 



k^to 



+k' 



9 hi J 



(B3) 



9C/ dU 



>^ ki • ka dU dU 



where we used the equation for the Jacobian PT|) and the equahty (|Bip . After that we have the equation: 



92 [/ 



1+ V ^1-^2 9t/ 



(B4) 



With help the relationship (|Bip we can write 



1 9 



2 5Pk, 



9 In J 1 kik2 9 In J 



= 0, 



(B5) 



because the expression in the square brackets coincides with the equation (ITT|) . We can see that the wave function 
(jBl|l satisfies to the Schrodinger equation (fT5|) . 

Let's obtain normalization constant i3. Let the wave function of a ground state ■0 (|23p is determined by the following 
way in RPA: 



ip = A exp 



1^ 1 

k#o 



(B6) 



where was taken into account. The constant A is obtained from the normalization condition ([7]) in the thermo- 
dynamical limit iV — ^ 00: 



' />oo />oo ^ / -1 



Let's write the Jacobian in RPA 

J = C exp 

The constant C is obtained from the condition ([5]) in the thermodynamical limit N ^ 00: 



(B7) 



k/0 



(B8) 



(B9) 



Then 



- 1 



(BIO) 
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If interaction between particles is absent, then S{k) = 1 and the wave function of a system ijj transforms to the wave 
function of a system of free bosons ^i^=o = 1/ VV^ , and the normahzation constant is 

B = 1/VV^. (Bll) 
The constant B is the same for higher powers of expansion of the wave function and Jacobian. 
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